A Time Delay for the Largest Gravitationally Lensed Quasar: 

SDSS J1004+4112 
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ABSTRACT 

We present 426 epochs of optical monitoring data spanning 1000 days from December 2003 
to June 2006 for the gravitationally lensed quasar SDSS J1004+4112. The time delay between 
the A and B images is A<ba = 38.4 ± 2.0 days (A% 2 = 4) in the expected sense that B leads A 
and the overall time ordering is C-B-A-D-E. The measured delay invalidates all published mod- 
els. The models failed because they neglected the perturbations from cluster member galaxies. 
Models including the galaxies can fit the data well, but strong conclusions about the cluster mass 
distribution should await the measurement of the longer, and less substructure sensitive, delays 
of the C and D images. For these images, a delay of A^cb — 681 ± 15 days is plausible but 
requires confirmation, while delays of A<cb > 560 days and A£ad > 800 days are required. We 
clearly detect microlensing of the A/B images, with the delay-corrected flux ratios changing from 
rriB — rriA = 0.44 ± 0.01 mag in the first season to 0.29 ± 0.01 mag in the second season and 
0.32 ± 0.01 mag in the third season. 
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1. Introduction 

The wide-separation lensed quasar SDSS J10044 
was discovered in the Sloan Digital Sky Survey 
search for lenses (Inada et al. 2003; Oguri et al. 
2004; Sharon et al. 2005; Wambsganss 2003). 
The lens consists of at least four images of a red- 
shift z s = 1.734 quasar whose ~ 15'.'0 Einstein 
ring diameter is created by a redshift zi — 0.68 
cluster. The cluster has been characterized with 
X-ray observations (Ota et al. 2006) and there are 
additional multiply imaged arcs formed from still 
higher redshift background galaxies (Sharon et al. 
2005). There is also strong evidence for a fifth, 
lensed image of the quasar located near the center 
of the brightest cluster galaxy (Inada et al. 2005), 
which in combination with a future velocity dis- 
persion measurement for the galaxy will strongly 
constrain the central mass distribution of the lens 
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(e.g. Sand et al. 2002, but see Dalai & Keeton 
2003). Thus, it is not only feasible to cleanly com- 
pare X-ray and lensing mass distributions in this 
galaxy cluster, but it may also be possible to test 
the cosmological model by measuring the increase 
of the Einstein radius with source redshift due to 
the Dls/Dos distance ratio scaling of the lens 
deflection (Soucail et al. 2004). 

That the source is a time- variable quasar offers 
further and unique opportunities for this cluster 
lens. First, the time delay between the quasar 
images can be measured as a constraint on the 
mass distribution. In theory, the time delays de- 
termine the mean surface density near the im- 
ages for which the delay is measured (Kochanek 
2002), so the mass sheet (k) degeneracy of most 
cluster lensing measurements can be broken un- 
der the assumption that the Hubble constant is 
well-determined by other means. Several theoreti- 
cal studies of the time delays in SDSS J1004+4112 
(Oguri et al. 2004; Williams & Saha 2004; Kawano 
& Oguri 2006) have explored their dependence 
on the mean mass profile of the cluster, finding 
a broad range of potential delays. As we shall 
see, all these models are incorrect in their details 
because they neglected cluster member galaxies 
whose deflection scales are larger than the posi- 
tional constraints on the quasar images used in 
the models (see the discussion in Kccton et al. 
2000 on the failure of similar models for the clus- 
ter lens Q0957+561 and the general discussion in 
Kochanek 2005). Nonetheless, all these models in- 
dicate that the delay between the A and B images 
is relatively short (weeks) and that its value should 
indicate the magnitude of the much longer (years) 
delays of the C and D images. 

The second unique property of the lens is that 
microlensing of the quasar accretion disk by any 
stars in the cluster halo or small satellites near 
the images can be used as an added probe of the 
structure of the cluster (see Wambsganss 2006). 
Because the cluster has a higher velocity disper- 
sion (700 km/s) than a typical galaxy lens (~ 
200 km/s), the microlensing time scales in this sys- 
tem may also be shorter than for a lens by about 
a factor of 3. There is already evidence for mi- 
crolensing from the time variability of the C IV 
1549 A line in image A that is not observed in im- 
age B (Richards et al. 2004, Lamer et al. 2006, 
Gomez- Alvarez et al. 2006), although recently 



Green (2006) has suggested that this could also 
be due to time variable absorption in the source 
quasar. 

For three years we have conducted an optical 
monitoring campaign to measure the optical vari- 
ability of this system. This has proved more chal- 
lenging than desired because the quasars are some- 
what faint for monitoring with available telescopes 
and modest exposure times. However, we have 
succeeded both in measuring the A/B time delay 
and clearly detecting microlensing of the optical 
continuum of the quasar. In §2 we present the 
data from the monitoring campaign for the four 
bright lensed quasar images. In §3 we determine 
the A/B time delay, discuss the presence of mi- 
crolensing in the system, and place constraints on 
the long delays between the close image pair A 
and B and the fainter images C and D. In §4 we 
discuss the failure of existing models for the sys- 
tem and introduce a simple successful model that 
includes the perturbations of cluster galaxies, and 
we conclude in §5. 

2. Data 

The photometric monitoring observations pre- 
sented here took place between December 2003 
and June 2006. The bulk of data were taken 
with the 1.2m telescope at Fred Lawrence Whip- 
ple Observatory on Mount Hopkins using the 
4Shootcr (R-band, 93 epochs, 0'.'66 pixels), Mini- 
cam (SDSS r-band, 74 epochs, 0'.'604 pixels), and 
Kcplercam (SDSS r-band, 91 epochs, / .'672 pixels, 
plus 4 epochs in R-band) during the first, second 
and third season, respectively. Additional data 
were obtained with the Apache Point Observa- 
tory (APO) 3.5m telescope using SPICam (SDSS 
r-band, 9 epochs, 0'.'282 pixels), the MDM 2.4m 
Hiltner telescope using the RETROCAM (Mor- 
gan et al. 2005, SDSS r-band, 27 epochs, 0"259 
pixels), 8K (R-band, 12 epochs, 0'/344 pixels), 
Tcmpleton (R-band, 8 epochs, 0'.'275 pixels) and 
Echelle (R-band, 3 epochs, 0'.'275 pixels) detec- 
tors, the MDM 1.3m McGraw-Hill telescope using 
the Templeton detector (R-band, 6 epochs, 0'.'508 
pixels), the Palomar Observatory 1.5m telescope 
using the SITe detector (R-band, 13 epochs, 0'.'379 
pixels), the Wise Observatory 1.0m telescope with 
the Tektronix (R-band, 30 epochs, 0'.'696 pixels) 
and TAVAS (clear, 53 epochs, 0'.'991 pixels) de- 
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Fig. 1. — The r-band image obtained with Mini- 
cam on January 16, 2005. The 3.3 arcmin x 
3.5 arcmin field shows the lensed images of SDSS 
J1004+4112 and the five reference stars SI, S2, 
S3, S4 and S5 used for the PSF. 



tectors, and the WIYN 3.5m telescope using the 
WTTM (SDSS r-band, 3 epochs, 0'/216 pixels) 
detector. The combined data set consists of 426 
epochs. 

In Figure 1 the quasar images are labeled A, B, 
C and D, following the notation by Inada ct al. 
(2003). The (non- variable) reference stars used 
for flux calibration and building the PSF are SI, 
S2, S3, S4 and S5. The small panels in Figure 
2 show snapshots of the four bright quasar im- 
ages at three different observing epochs, in March 
2004, May 2005 and March 2006. These images 
illustrate how images A and B slowly faded dur- 
ing the course of the three seasons, while image D 
became significantly brighter. The galaxies of the 
lensing cluster are not detectable in the individual 
observations, except for the bright galaxy close to 
image D (Gl in Oguri et al. 2004). The candidate 
fifth quasar image, E, lies near the center of this 
galaxy (Inada et al. 2005). 

The data were fitted using the methods of 
Kochanek et al. (2006) for HE 0435-1223. Re- 
gions around each of the quasar images and the 



Fig. 2.— The panels show 23" x 23" insets on the 
four bright quasar images at three different epochs 
separated by about one year (for nomenclature cf. 
Fig. 1). The faint image in their middle is the 
bright galaxy belonging to the lensing cluster. 

"standard" S1-S5 stars (see Fig. 1) are fitted to 
determine the relative fluxes and the structure 
of the PSF. For each filter, the star SI was de- 
fined to have unit flux while the fluxes of the re- 
maining stars S2, S3, S4 and S5 were adjusted to 
this calibration standard based on all the available 
epochs of data for each filter. The relative fluxes 
of the standard stars depend on the filter, with ra- 
tios of 1.0:0.439:0.360:0.130:0.0583 for the R-band, 
1.0:0.334:0.329:0.0937:0.0613 for the SDSS r-band, 
and 1.0:0.63:0.64:0.39:0.20 for the clear filter. In 
the WIYN/WTTM, MDM 2.4m/8K and MDM 
2.4m/Templcton data, the star SI frequently is 
too close to saturation for use, so its weight in 
the fits is greatly reduced. It was not necessary 
to further subdivide the calibrations for the indi- 
vidual detectors given the overall quality of the 
photometry, as the average calibration offsets be- 
tween detectors were well under 0.01 mag. We 
then matched the R-band and clear observations 
to the r-band observations using the quasar light 
curves themselves. For each R/clear epoch brack- 
eted by r-band observations within 1 week, we in- 
terpolated the r-band observation to the epoch of 
the other band and computed the mean offset be- 
tween the light curves. Offsets of 0.043±0.006 mag 
and 0.250 ± 0.011 mag must be added to the R- 
band and clear magnitudes respectively to match 
them to the r-band data. 

Figure 3 shows the resulting light curves for im- 
ages A-D over the three observing seasons and 
Table 1 presents the photometry for images A and 
B. They span a time period of 1000 days from 
December 2003 to June 2006 with two seasonal 
gaps of approximately 100 days during the period 
from July to October. SDSS J1004+4112 is a rel- 
atively faint quasar for monitoring with lm-class 
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Fig. 3.— Light curves of the A-D images of SDSS J1004+4112 from December 2003 to June 2006. Images C 
and D have been offset by 0.3 and 0.6 mag because they would otherwise overlap with each other and with 
image B in the third season. 
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telescopes, and the image quality of the FLWO 
and WISE telescopes is poor. As a result, the 
noise in many of the measurements is relatively 
large compared to the variability amplitude. On 
the other hand, our sampling cadence is quite high, 
so the overall statistical power of the data is very 
good, with a mean sampling rate of once every two 
days while the source is visible. All four images 
vary by about 0.5 mag, with the more than 1 mag 
brightening of image D being the largest change 
during the three seasons. For the purposes of mea- 
suring the A/B time delay, the most interesting 
features are the minima in the B light curve near 
days 3150 and 3750 in the first and third seasons 
respectively, and the corresponding features in the 
A light curve roughly 40 days later. The second 
season shows no obvious features that can be used 
to measure the delay. The second important point 
to note is that the A/B flux ratio has changed 
significantly between the first and third seasons, 
indicating that microlensing is occurring in this 
system as has been previously suggested by vari- 
ations in the C IV emission line profile (Richards 
et al. 2004, Lamer et al. 2006, Gomez- Alvarez et 
al. 2006). 

3. The Time Delay 

Model predictions for the time delay of the close 
image pair A and B are a few weeks (Oguri et al. 
2004; Williams & Saha 2004, Kawano & Oguri 
2006) and therefore should be measurable within 
each season of the light curves. Of the many 
techniques for calculating time delays from light 
curves (e.g. Gil-Merino et al. 2002, Pelt et al. 
1994, Press, Rybicki & Hewitt 1992, Kochanek et 
al. 2006), we will apply three. The three meth- 
ods produce mutually consistent results, but we 
will adopt the Kochanek et al. (2006) polynomial 
method for our standard result because it natu- 
rally includes the effects of microlensing on the 
delay estimate. As is clear from the light curves, 
image B leads image A, so the delay ordering of 
the images is C-B-A-D-E. We conclude with a dis- 
cussion of the longer C and D image time delays. 

For our analysis of the A/B delay we treated 
the data in Table 1 as follows. If the goodness of 
fit of the photometric model to an image had a 
X 2 statistic larger than the number of degrees of 
freedom Ndof (see Table 1), we rescaled the pho- 



tometric errors for that image by (x 2 /Ndof) on 
the grounds that having x 2 > ^dof meant that the 
uncertainties were underestimated. For the time 
delay estimates we dropped the 16 points marked 
in Table 1 that were more than 3cr from the best 
fitting models. We also repeated the time delay 
estimates excluding all points with rescaled pho- 
tometric errors larger than 0.1 mag, finding no sig- 
nificant changes. 

3.1. Simple x 2 minimization 

The simplest approach to the delay measure- 
ment problem is to take the observed light curves 
A(U) and B(ti) and cross-correlate them with lin- 
early interpolated light curves a{t) and b{t) for the 
other image. We assume that that the light curves 
of the two images are the same except for a time 
delay r and a magnitude offset m(r). In prac- 
tice, we use a different magnitude offset for each 
season to partially compensate for the effects of 
microlensing. Based on this assumption we can 
calculate the time delay by minimizing the devia- 
tions from m(r) for each pair (A(ti), b(ti — r)) and 
(a(ti + t), B{ti)) by a fit statistic 

X 2 ( T ) _ _J (A(U)~b(t t -T)+m(r)) 2 

N dof (T) 27V do/ (r) ^+< t 

N(t) 

that is symmetric as to which image is being in- 
terpolated. The errors in the observed magnitudes 
are a as and as.i and the errors in the interpolated 
magnitudes are a a> t and <7b ;t . The fit is carried out 
only where the light curves overlap (i.e. excluding 
the season gaps), so the number of data points 
used N(t) depends on the delay r. 

Fig. 4 shows the results for the three seasons 
separately and for the combined light curve. An- 
alyzed separately, the first and third seasons show 
minima at 31 and 40 days respectively, while there 
is no clear minimum for the second season due to 
the lack of significant features in the light curve. 
For the joint analysis of all three seasons we al- 
lowed for an independent value of m(r) within 
each season to model the changes in the flux ratios 
due to microlensing. The analysis of the combined 
data yields a delay of 41±5 days. 
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Fig. 4. — Results of the x 2 minimization between 
the two time series A and B for the first (red dot- 
ted line), second (green dashed line), third (blue 
dash-dotted line) and combined observing seasons 
(black solid line). 

3.2. The Dispersion Method 

One potential weakness of the simple \ 2 
method is the need for interpolation. As our 
second approach we apply the dispersion spec- 
tra method developed by Pelt et al. (1994, 1996) 
to avoid the interpolation. Instead, a combined 
light curve C (t) is constructed by shifting the data 
points of one image in magnitude (m(r)) and time 
(t) and combining them with the data points of 
the other image 



C(t fc ) 



Aj, if t k — ti 

Bj — to(t), if t k — tj + r 



(2) 



where k = 1, .., N and N = Na + N B - The time 
delay r is estimated by minimizing the dispersion 
spectrum 



D 2 (r) 



N-l 

E S k W k G k (C k+1 -C fe ) 2 
™R — — ( 3 ) 

m(r) 



2 J2S k W k G k 

k=l 



where the W k = {a\ + cr^ +1 ) _1 are the statisti- 
cal weights of the data, G k = 1 if the points k 
and k + 1 come from different images (A/B) and 
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Fig. 5. — Dispersion spectra for the first (red dot- 
ted line), second (green dashed line), third (blue 
dash-dotted line) and combined observing seasons 
(black solid line). 

G fe = otherwise (A/A or B/B), and S k = 1 if 
\t k+ i — tk\ < S and S k = otherwise. We use 
a decorrelation time scale of S = 3 days, but our 
results depend little on the exact choice. The re- 
sults are shown in Figure 5 for both the individual 
seasons and the combined data. We again used 
independent estimates of Am for each observing 
season to compensate for the effects of microlens- 
ing. We find 34 and 39 days for the first and third 
seasons, 39 days for the combined data, and no 
significant minimum using only the data from the 
second season. 

We estimated the errors using the resampling 
procedure of Pelt et al. (1994). The combined 
light curve C k was smoothed for each time de- 
lay using a 7-point median filter surrounding each 
point. Residuals relative to the original data 
were then reshuffled randomly to create artificially 
noisy combined light curves. Time delays for a set 
of 1000 such light curves were determined by cal- 
culating the dispersion spectra, leading to the dis- 
tribution of minimum dispersion estimates shown 
in Fig. 6. If we define the uncertainties by the 
range about the median encompassing 68% of the 
random trials, we estimate that the uncertainty in 
the time delay is ±6 days. 

3.3. The Polynomial Method 

The clear indication of microlcnsing effects 
means that corrections for microlensing are re- 
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Fig. 6. — Results of the resampling procedure in 
the dispersion method. From the width of the dis- 
tribution we estimate the uncertainty for the time 
delay measurement in the dispersion spectra. 

quired to determine an accurate time delay. Both 
the x 2 an d minimum dispersion methods treated 
the flux ratios between the images within each 
season as a constant. Either method could be 
modified to allow for more complex microlensing 
variations, but for our final analysis we will use 
the polynomial fitting method of Kochanek et al. 
(2006) since it can most easily incorporate the ef- 
fects of microlensing on both the delays and their 
uncertainties. 

In the Kochanek et al. (2006) polynomial 
method, the time variations of the source are mod- 
eled as a Legendre polynomial of order iV src , and 
the time variations due to microlensing are mod- 
eled as a Legendre polynomial of order in each 
of the three seasons. The amplitudes of the coef- 
ficients of the source polynomial are weakly con- 
strained to match the structure function measured 
for SDSS quasars by Vanden Berk et al. (2004). 
The polynomial orders are determined by using 
the F-test to indicate which polynomial order no 
longer leads to statistically significant improve- 
ments in the fits. We used polynomial orders of 
N src = 20, 40, and 60 and = 0, 1 and 2. 
The microlensing polynomial orders correspond to 
using a constant flux ratio, a linear trend or a 
quadratic trend for each season. Based on the 
F-test, the improvement in the fit to the data is 
significant when jumping from N src = 20 to 40 
and from = to 1 (from constant flux ratios 



in each season to linear trends), but not for any 
of the higher-order models. The delays for all the 
cases are consistent with each other given their 
uncertainties, so we will adopt the result for the 
N src = 60, N^ = 3 model, At B A = 38.4± 1.0 days 
(Ax 2 = 1, ±2.0 days at A\ 2 = 4). Using higher 
than necessary polynomial orders should be con- 
servative and overestimate the uncertainties in the 
time delay. The overall fit has % 2 = 718 for 
Ndof = 663 degrees of freedom. 

In this model, the mean magnitude differences 
between A and B for the three seasons are 0.439 ± 
0.008, 0.292 ± 0.012 and 0.321 ± 0.008 mag, with 
seasonal gradients of -0.10 ± 0.03, -0.27 ± 0.04 
and —0.11 ±0.03 mag/year and second derivatives 
of -1.2 ± 0.4, 0.6 ± 0.5 and -0.7 ± 0.3 mag/year 2 
respectively. Thus, microlensing is clearly present, 
as expected from the visible structure of the A and 
B light curves. The need to model the microlens- 
ing as more than a seasonal change in the flux ra- 
tio means that the polynomial models fit the data 
considerably better than the first two methods, 
which is one reason for the significantly smaller 
formal uncertainties in the delay. Using only the 
higher precision data points has a negligible effect 
on the delays or the inferred level of microlens- 
ing. Fig. 7 shows the estimated source light curve 
as compared to the data, and Fig. 8 shows the 
inferred level of microlensing variability. We can 
only measure the differential microlensing between 
A and B, and the choice of assigning it to image B 
is an arbitrary one which does not affect the time 
delay estimate. 

3.4. Constraints on the long delays 

The model predictions for the long time delays 
between the close image pair A and B and the 
fainter images C and D are very uncertain. For 
example, Oguri et al. (2004) found an approxi- 
mate scaling relation of AicD/AteA = 143 ± 16 
for their models, which would imply a 15 year C/D 
time delay given our results for the A/B time de- 
lay. On the other hand, Williams & Saha (2004) 
found delay estimates of order AicB ~ 400 days 
and At ad ~ 600 days, albeit with a large scatter 
(about ±200 days). As we will discuss in §4, these 
studies use simplified mass models that are of only 
limited use for calibrating our expectations. 

Empirically, with our 1000 day time span for 
the light curves we can test for image C and D de- 
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Fig. 7. — The overlapping A (squares) and B 
(triangles) light curves after shifting by the time 
delay and subtracting the estimated microlensing 
variability. The curve is the best fit N src = 60 
source light curve model. There are three inde- 
pendent fits, one for each season, connected by a 
line through the seasonal gaps. The stars are the 
points that were masked in the time delay analysis. 
To minimize the confusion we are showing only the 
data with uncertainties smaller than 0.1 mag and 
include examples of their mean error bars in the 
lower left corner. 

lays of ±1000 days. We did so by matching the C 
or D light curve to the combined A/B light curves 
using the polynomial method. Since the overall 
behavior of the A and B light curves during the 
first and second season is mainly decreasing or flat 
while the light curve of image C shows an increase 
in the first season (Fig. 3), the time delay between 
C and B (with C leading) must be larger than 
560 days. Assuming C is leading, there is a min- 
imum near 681 ± 15 days which corresponds to 
aligning the minimum observed in the first sea- 
son for C with that observed in the last season for 
A/B. Due to the very different shapes of the light 
curves of images A/B and D there is no obvious 
solution over the whole observed time span. The 
only possibility would be to match the plateau in 
the third-year data of image D with the initial 
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Fig. 8. — The inferred microlensing variability. 
We assigned the microlensing variability to image 
B, leaving the model for image A as a constant. 
This is an arbitrary choice that has no effect on the 
time delay - all we can really measure is the differ- 
ential microlensing between the two images. The 
curves show the constant model used for A, the 
quadratic model (N^ = 3) for B and only the data 
with uncertainties smaller than 0.1 mag. There 
are three independent fits, one for each season, 
connected by a line through the seasonal gaps. 

portions of the first-year data of images A/B, but 
there is no good candidate minimum in the good- 
ness of fit. Therefore, we conclude that the time 
delay between A and D is larger than 800 days. 

4. Models and Interpretation 

We modeled the system using lensmodel (Kee- 
ton 2001) and the same component positions as 
were used by Kawanao & Oguri (2006) and In- 
ada et al. (2006). We fitted all five quasar im- 
ages assuming astrometric minimum uncertainties 
of 0'.'003 and 20% flux uncertainties that are com- 
parable to the observed amount of microlensing. 
We used the accurate but slow image plane fit- 
ting method, and the Hubble constant was fixed 
at H = 72 km s" 1 Mpc -1 . The brightest clus- 
ter galaxy was modeled as an ellipsoidal de Vau- 



couleurs model with a major axis effective radius 
of R e = 4'.'5, an ellipticity of 0.18±02 and a major 
axis position angle of— 12° ±5° based on fits to 
the CASTLES project's Hubble Space Telescope 
(HST) NICMOS H-band image of the system. The 
cluster halo was modeled as an ellipsoidal NFW 
model with a break radius of r s = 40'.'0 based on 
the mass model for the X-ray emission by Ota et 
al. (2006). We assumed priors on the ellipticity 
of the halo of 0.18 ± 0.05, no prior on its major 
axis position angle, and a prior on the external 
shear of 7 = 0.05 ± 0.05. We also imposed hard 
limits on the ellipticity and position angle of the 
central galaxy (0.15 to 0.25 and -22° to -2°), the 
galaxy position (±0'.'3 in each coordinate), the el- 
lipticity of the halo (0.0 to 0.5), the position of the 
halo (±3'.'0 relative to the central galaxy) and the 
shear (0 < 7 < 0.25). 

We first ran a model sequence based on simply 
adding a halo to the central galaxy. We started 
by fitting a de Vaucouleurs model with no halo to 
get the mass scale needed for the central galaxy in 
the absence of a halo. Then we fitted a series of 
models with the mass of the central galaxy fixed 
to a fraction < / < 1 of its value in the no halo 
model. We ran the series both with and without 
the putative fifth quasar image. In general, the re- 
sults are poor. The best fits to the image positions 
are obtained for / ~ 0.1. The time delays for these 
models strongly disagree with our measurement in 
the sense that the model A-B delays are too short 
(<~ 15 days for / = 0.1). Producing a longer delay 
requires a model with a lower surface density near 
the images, since the time delays of these simple 
models are roughly proportional to 1 — (k) where 
(k) is the mean surface density in the annulus be- 
tween the images (see Kochanek 2002). However, 
the models with / ~ 0.7 and a low surface density 
which fit the delay correctly, fit the images poorly 
and have cllipticities for both the galaxy and the 
halo that are driven to their maximum permitted 
values because a side effect of lowering the sur- 
face density is to increase the required ellipticity 
(see Kochanek 2005). That these simple models 
fit our delay measurement poorly is not surprising 
since the published results based on these simple 
model classes 1 never produced a delay as long as 

lr The non-parametric models of Williams & Saha (2004) arc 
roughly comparable in their overall structures. We note in 
passing that the discrimination between radial mass pro- 



our measured value. 

The fundamental problem with this model, and 
all the preceding models of Oguri et al. (2004), 
Williams & Saha (2004), and Kawano & Oguri 
(2006), is that they neglect or poorly represent 
the substructure in the potential due to the pres- 
ence of the other cluster galaxies. Many of these 
galaxies have deflection scales that are enormous 
compared to the astrometric uncertainties in the 
image positions, and as we painfully learned over 
20 years of modeling Q0957 +561, astrometric un- 
certainties can be imposed to no greater accuracy 
than the deflection scales of the most massive ne- 
glected components of the mass distribution (see 
Keeton et al. 2000, Kochanek 2005). In short, 
the model sequence we just considered, as well as 
all published models of this system, was virtually 
guaranteed to be quantitatively incorrect. 

At a minimum, the model needs to include 
galaxies whose deflections cannot be trivially mim- 
icked by rescaling the mass of the central galaxy 
and modifying the external shear. We used Sex- 
tractor (Bertin & Arnouts 1996) to determine the 
positions and fluxes of the galaxies in the CAS- 
TLES HST/ ACS I-band image of the cluster. We 
assumed the galaxies had critical radii that scaled 
with the square root of their flux (i.e. SIS mod- 
els obeying a Tully-Fishcr relation), and added 
the 11 most important galaxies within 20'.'0 of the 
main lens galaxy as circular pseudo-Jaffe mod- 
els (p oc r~ 2 [r 2 + a 2 ] -1 ) with a break radius of 
a = l'.'O. We required that they have mass scales 
(Einstein radii) in the range 0705 < b < 2'.'0 and 
kept their positions fixed. 2 We did not attempt 
to force a correlation between flux and Einstein 
radius as the scatter in the relation is fairly large 
(Rusin et al. 2003). Figure 9 shows the posi- 
tions of these galaxies relative to the image po- 
sitions and the cluster center. We then ran the 
same sequence of models for the central galaxy 

files observed in the non-parametric models is purely an 
artifact of the priors used in the analysis - there is a math- 
ematical degeneracy that makes it impossible to use the 
positions of images A-D to determine the radial mass pro- 
file without the further assumptions supplied by the priors 
(see Kochanek 2005). In our case, adding image E partly 
breaks the degeneracy. 
2 In order of increasing RA, they were the galaxies located 
at (-14'.'8,-5'.'5), (-12"2,3'.'7), (-12'.'2,13'.'5), (-9"2,2'.'5), 
(-9'.'5,18"6), (-6'.'5,11'.'3), (-4V4,-0^'6), (-2'.'8,-0'.'2), 
(-2'.'1,11{'9), (-2('8,14'.'l) and (l('4,0'/5) from image A. 
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Fig. 9. — A lens model generally consistent with 
the data. The open squares show the positions 
of the 5 lensed images A-E. The solid lines show 
the critical curves of the lens model, whose con- 
siderable, but realistic, complexity is due to the 
presence of the main lens galaxy (filled square), an 
offset halo (filled triangle), and 11 smaller galaxies 
(large crosses) . The model fits the time delays and 
the image positions nearly perfectly. 

and halo. These models have no difficulty fitting 
both the time delay and the image positions with 
reasonable parameters and a dark-matter domi- 
nated cluster model (/ ~ 0.1). Given the sensitiv- 
ity of the A/B delay to substructure, it is probably 
premature to use the time delays as a strong con- 
straint on the structure of the cluster. Reasonable 
models predict B/C delays of order 450 to 1000 
days, suggesting that the roughly 680 day solution 
in §3.4 may well be correct, and that the A/D de- 
lays are of order 5-7 years. These longer delays 
should be much less sensitive to the perturbations 
from galaxies and will provide a better basis for 
studying the cluster. 

Figure 9 shows the critical line structure of an 
illustrative model. Since the best models lead to 
B/C delays quite close to the value of 681 days 
found by matching the minimum in the first sea- 
son for C with those in the last season for A/B, we 
added it as a constraint. The model fits the im- 



age positions and the two time delays well but not 
perfectly, and there is a modest (x 2 = 11.5) ff.'l 
westward offset of the lens galaxy suggesting that 
there are some remaining issues for fitting image 
E. In this model, the D image time delay relative 
to A and B is approximately 5.7 years and the cen- 
tral galaxy has 10% of the mass it would require 
in the absence of the dark matter halo. The Ein- 
stein radii of the smaller galaxies range from 0'.'05 
to l'/4, which are plausible mass scales. The offset 
between the main lens galaxy and the dark mat- 
ter halo found in earlier models still seems to be 
required. It is not our present intent to conduct a 
full model survey, but to emphasize the need for 
realistic models. 

5. Summary and Conclusions 

We have measured the A/B time delay of SDSS 
J1004+4112 to be (38.4 ± 2.0) days {A X 2 = 4), 
which fixes the overall time ordering of the images 
to be C-B-A-D-E. While this is the time ordering 
predicted in published models (Oguri et al. 2004, 
Williams & Saha 2004, Kawano & Oguri 2006) it 
is significantly longer than the delays predicted by 
these models. The cause of the discrepancy is that 
the previously published models overly simplified 
the mass distribution by neglecting the deflections 
generated by the cluster member galaxies. 

Models including the eleven most important 
galaxies can simultaneously fit the A-E image po- 
sitions and the measured A/B time delay with rea- 
sonable parameter values. Modelers of this sys- 
tem need to remember the lesson of Q0957+561 
- model constraints that are applied more tightly 
than the deflection scale of the most massive, ne- 
glected components of the lens lead to incorrect 
results (Keeton et al. 2000). We note that Sharon 
et al. (2005) also needed to include some of the 
member galaxies in order to model the higher red- 
shift lensed arcs, but made no predictions for the 
time delay. 

Fortunately, the A/B delay should be the most 
sensitive of the delays to the effects of cluster 
galaxies because it is a merging image pair. The 
longer delays for the C and D images relative to A 
and B should be less affected by substructure, so 
their measurement should provide constraints on 
the cluster halo properties that are less sensitive to 
the member galaxies. At present we cannot claim 
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a measurement of these longer delays. A lower 
bound on the delay A£ad > 800 days is consistent 
with our models, which predict delays of 5-7 years 
for this image pair. The shorter C/B delay is at 
least AicB > 560 days but there is a possible delay 
of AtcB — 680 ± 15 days that should be confirmed 
or rejected during the next observing seasons and 
is consistent with our models. 

We have also clearly detected microlensing vari- 
ability in the A/B images, with changes of order 
0.15 mag in the A/B flux ratio over the course 
of the three observing seasons. This result pro- 
vides strong evidence that the differential changes 
in the A/B emission line profiles are also due to 
microlensing (Richards ct al. 2004, Lamer et al. 
2006, Gomez- Alvarez et al. 2006) rather than 
variable absorption in the source (Green 2006). 
The microlensing time scales in SDSS J1004+4112 
should be relatively shorter than in most single 
galaxy lenses because the internal velocities of the 
cluster are about 3 times higher than those of a 
galaxy. While the flux ratio changes in the optical 
continuum are modest, we would expect to find 
significantly larger effects at shorter wavelengths 
where the source size should be more compact. 
There is already some evidence for this from the 
X-ray flux ratios measured by Ota et al. (2006) 
and Lamer et al. (2006). A campaign to monitor 
this system in X-rays would both allow us to study 
the size of the X-ray emission region and provide 
the added data on the emission from the cluster 
needed to provide a precision comparison of the 
mass distributions estimated using X-ray data and 
lens models. Such careful tests will be essential if 
measurements of the increase of the Einstein ra- 
dius of the cluster with source redshift based on 
the surrounding multiply imaged arcs are to be 
used as a new test of the cosmological model as 
proposed by Soucail et al. (2004) and Sharon ct 
al. (2005). 
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Table 1 

Light Curves for SDSS J1004+4112 A&B 



HJD x 2 /Ndof Image A Image B Observatory Detector 
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Table 1 — Continued 



HJD X 2 /Ndof Image A Image B Observatory Detector 
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3 


.563 


± 


0. 


.052 


Wise 


m A t r A i~i 

TAVAS 


3064 


.468 


35.89 


3 


.823 


± 


0. 


.054 


5 


.480 


± 


0. 


.226 


Wise 


Tektronix 


3064 


.773 


1.01 


3 


.112 


± 





,017 


3 


.603 


± 


0. 


.026 


FLWO 


4Shooter 


3064 


.884 


1.41 


3 


.121 


± 


0. 


.007 


3 


.570 


± 


0. 


.008 


"TV T T "\ rv I 

MDM 


Echelle 


3065 


.480 


12.90 


3 


,191 


± 





.057 


3 


.958 


± 


0. 


,112 


Wise 


Tektronix 


3065 


.805 


2.91 


3 


.099 


± 


0. 


.008 


3 


.547 


± 


0, 


Oil 


FLWO 


4Shooter 


3073 


.862 


0.59 


3 


.086 


± 


0. 


.023 


3 


.547 


± 


0, 


,034 


FLWO 


4Shooter 


3075 


.464 


1.03 


3 


.091 


± 


0. 


.047 


3 


.740 


± 


0. 


,083 


Wise 


Tektronix 


3078 


.837 


1.12 


3 


.077 


± 


0. 


.012 


3 


.584 


± 


0, 


,018 


FLWO 


4Shooter 


3079 


.832 


1.77 


3 


.048 


± 


0. 


.013 


3 


.531 


± 


0, 


,019 


FLWO 


4Shooter 


3079 


.848 


11.73 


3 


.026 


± 





.007 


(3 


.471 


± 


0. 


.007) 


TV /T~T~V TV IT 

MDM 


8K 


3080 


.822 


1.01 


3 


.077 


± 


0. 


,013 


3 


.566 


± 


0. 


.019 


FLWO 


4Shooter 


3081 


.327 


2.06 


3 


.005 


± 


0. 


.062 


3 


.454 


± 


0. 


.099 


Wise 


TAVAS 


3081 


.789 


3.07 


3 


.107 


± 


0. 


,011 


3 


.587 


± 


0. 


.015 


FLWO 


4Shooter 


3082 


.785 


1.69 


3 


.072 


± 


0. 


,011 


3 


.576 


± 


0. 


,016 


FLWO 


4Shooter 


3083 


.934 


1.01 


3 


.113 


± 


0. 


.022 


3 


.517 


± 


0. 


,031 


FLWO 


4Shooter 


3085 


.793 


1.45 


3 


.087 


± 


0. 


.011 


3 


.539 


± 


0. 


,015 


FLWO 


4Shooter 


3087 


.662 


3.09 


3 


.012 


± 


0. 


.061 


3 


.498 


± 


0, 


,061 


WIYN 


tit m m Ti ,t 

WTTM 


3087 


.707 


2.40 


3 


.077 


± 


0. 


.011 


3 


.534 


± 


0, 


,016 


FLWO 


4Shooter 


3088 


.734 


2.01 


3 


.063 


± 


0. 


.011 


3 


.563 


± 


0. 


,016 


FLWO 


4Shooter 


3090 


.735 


0.63 


3 


.088 


± 


0. 


.024 


3 


.587 


± 


0, 


,036 


FLWO 


4Shooter 


3091 


.772 


0.97 


3 


.080 


± 


0. 


,014 


3 


.534 


± 


0. 


.020 


FLWO 


4Shooter 


3092 


.821 


0.65 


3 


.069 


± 





,017 


3 


.527 


± 


0. 


.025 


FLWO 


4Shooter 


3093 


.790 


0.72 


3 


.068 


± 





,014 


3 


.546 


± 


0. 


,021 


FLWO 


4Shooter 


3094 


.718 


0.80 


3 


.022 


± 


0. 


.020 


3 


.510 


± 


0. 


.031 


FLWO 


4Shooter 


3095 


.819 


0.67 


3 


.080 


± 


0. 


.024 


3 


.516 


± 


0. 


,035 


FLWO 


4Shooter 


3100 


.412 


0.36 


3 


.085 


± 


0. 


,093 


3 


.529 


± 


0. 


.137 


Wise 


Tektronix 


3101 


.657 


0.42 


3 


.021 


± 


0. 


.026 


3 


.461 


± 


0, 


,038 


FLWO 


4Shooter 


3104 


.747 


0.99 


3 


.066 


± 


0. 


,017 


3 


.534 


± 


0. 


,025 


FLWO 


4Shooter 


3107 


.686 


1.31 


3 


.027 


± 


0. 


.755 


3 


.486 


± 


0, 


,755 


WIYN 


WTTM 


3107 


.694 


0.66 


3 


.065 


± 


0. 


.015 


3 


.591 


± 


0, 


,024 


FLWO 


4Shooter 


3108 


.285 


1.84 


3 


.050 


± 


0. 


.031 


3 


.606 


± 


0, 


,049 


Wise 


Tektronix 


3108 


.729 


1.64 


3 


.089 


± 


0. 


,012 


3 


.544 


± 


0. 


,016 


FLWO 


4Shooter 


3111 


.688 


1.61 


3 


.083 


± 


0. 


,013 


3 


.546 


± 


0. 


,018 


FLWO 


4 Shooter 


3113 


.335 


0.89 


3 


.071 


± 


0. 


,051 


3 


.566 


± 


0. 


,080 


Wise 


Tektronix 


3116 


.725 


1.37 


3 


.115 


± 





,013 


3 


.542 


± 


0. 


,017 


FLWO 


4Shooter 


3117 


.746 


1.36 


3 


.091 


± 


0. 


,014 


3 


.567 


± 


0. 


,020 


FLWO 


4Shooter 


3118 


.764 


0.82 


3 


.135 


± 


0. 


,014 


3 


.552 


± 


0. 


,019 


FLWO 


4Shooter 


3119 


.731 


0.76 


3 


.077 


± 


0. 


.014 


3 


.551 


± 


0, 


,020 


FLWO 


4Shooter 


3120 


.730 


0.70 


3 


.102 


± 


0. 


.015 


3 


.561 


± 


0, 


,022 


FLWO 


4Shooter 
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3122. 


.786 


0.55 


3 


.096 


± 





.028 


3 


.583 


± 





.043 


FLWO 


4Shooter 


3124. 


.718 


0.78 


3 


.083 


± 





.027 


3 


.614 


± 





.044 


FLWO 


4Shooter 


3125. 


.737 


0.52 


3 


.077 


± 





.034 


3 


.523 


± 





.051 


FLWO 


4Shooter 


3126 


.717 


0.58 


3 


.117 


± 





.027 


3 


.534 


± 





.039 


FLWO 


4Shooter 


3129 


.240 


0.72 


2 


.991 


± 





.057 


3 


.577 


± 





.094 


Wise 


Tektronix 


3129. 


.812 


0.58 


3 


.120 


± 





.035 


3 


.539 


± 





.051 


FLWO 


4Shooter 


3132 


.749 


2.09 


3 


.106 


± 





.008 


3 


.596 


± 





,012 


FLWO 


4Shooter 


3135. 


.717 


7.04 


3 


.016 


± 





.007 


3 


.544 


± 





.009 


MDM 


8K 


3136 


.659 


1.65 


3 


.077 


± 





.013 


3 


.625 


± 





.019 


FLWO 


4Shooter 


3137 


.719 


1.15 


3 


.045 


± 





.013 


3 


.627 


± 





.021 


FLWO 


4Shooter 


3140. 


.640 


1.58 


3 


.084 


± 





.013 


3 


.636 


± 





.019 


FLWO 


4 Shooter 


3143. 


.790 


0.63 


3 


,162 


± 





.029 


3 


.712 


± 





.047 


FLWO 


4Shooter 


3144. 


.677 


1.61 


3 


.089 


± 





,012 


3 


,659 


± 





.019 


FLWO 


4Shooter 


3145. 


.711 


1.55 


3 


.087 


± 





.014 


3 


.663 


± 





.023 


FLWO 


4Shooter 


3146. 


.745 


0.84 


3 


.093 


± 





,013 


3 


,636 


± 





.020 


FLWO 


4Shooter 


3147. 


.648 


1.42 


3 


.098 


± 





.015 


3 


,648 


± 





.023 


FLWO 


4Shooter 


3148. 


.653 


1.34 


3 


.078 


± 





.015 


3 


,684 


± 





.025 


FLWO 


4Shooter 


3149. 


.733 


1.37 


3 


.094 


± 





.013 


3 


.631 


± 





.020 


FLWO 


4Snooter 


3153. 


.660 


1.40 


3 


.060 


± 





.061 


3 


.631 


± 





.061 


WIYN 


WTTM 


3153. 


.692 


0.86 


3 


.120 


± 





.019 


3 


.661 


± 





,030 


FLWO 


4Shooter 


3154. 


.272 


0.33 


3 


.071 


± 





.059 


3 


.578 


± 





,093 


Wise 


Tektronix 


3154. 


.673 


0.66 


3 


.125 


± 





.023 


3 


.695 


± 





.037 


FLWO 


4Shooter 


3155. 


.673 


0.62 


3 


.106 


± 





.024 


3 


.613 


± 





.037 


FLWO 


4Shooter 


3156. 


.646 


0.72 


3 


.113 


± 





.020 


3 


.625 


± 





.031 


FLWO 


4Shooter 


3157. 


.670 


0.67 


3 


.100 


± 





.021 


3 


.688 


± 





.034 


FLWO 


4Shooter 


3158. 


.660 


1.04 


3 


.105 


± 





.017 


3 


,628 


± 





.026 


FLWO 


4Shooter 


3159. 


.643 


0.57 


3 


.145 


± 





.029 


3 


.679 


± 





.046 


FLWO 


4Shooter 


3161 


.736 


2.88 


3 


.116 


± 





.016 


3 


.658 


± 





.024 


MDM 


Templcton 


3162 


.658 


2.80 


3 


.150 


± 





.008 


3 


.653 


± 





,010 


MDM 


Templeton 


3163 


.650 


4.10 


3 


.159 


± 





.010 


3 


.653 


± 





,013 


MDM 


Templeton 


3164 


.651 


28.62 


3 


.168 


± 





.010 


3 


.646 


± 





,014 


MDM 


Templeton 


3165 


.656 


11.23 


3 


.171 


± 





.008 


3 


.623 


± 





.011 


MDM 


Templeton 


3166 


.650 


3.74 


3 


.142 


± 





.009 


3 


.615 


± 





.012 


MDM 


Templeton 


3167 


.665 


8.23 


3 


.136 


± 





.006 


3 


.621 


± 





.006 


MDM 


Templcton 


3168 


.663 


17.04 


3 


.152 


± 





.006 


3 


.636 


± 





.007 


MDM 


Templcton 


3169 


.662 


6.01 


3 


.171 


± 





.006 


3 


.635 


± 





.007 


MDM 


Templeton 


3169 


.674 


1.08 


3 


.164 


± 





.015 


3 


.642 


± 





.021 


FLWO 


4Shooter 


3170. 


.259 


0.28 


3 


.106 


± 





.047 


3 


.593 


± 





.073 


Wise 


Tektronix 


3170. 


.667 


0.71 


3 


.169 


± 





.011 


3 


.630 


± 





.016 


MDM 


Templcton 


3170. 


.673 


0.60 


3 


.174 


± 





.025 


3 


.687 


± 





.039 


FLWO 


4Shooter 


3171. 


.271 


0.33 


3 


.126 


± 





.048 


3 


.562 


± 





.072 


Wise 


Tektronix 
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3171. 


.661 


1 


.78 


3 


.185 


± 





.007 


3 


.634 


± 





.009 


MDM 


Tcmplcton 


3171. 


.666 


1 


.22 


3 


.119 


± 





.019 


3 


.635 


± 





.030 


FLWO 


4Shooter 


3172. 


.664 


1 


.68 


3 


.190 


± 





.007 


3 


.634 


± 





.008 


TV IT 1 "\ TV ,T 

MDM 


Tcmplcton 


3172. 


.668 





.73 


3 


.151 


± 





.016 


3 


.589 


± 





.022 


FLWO 


4Shooter 


3173. 


.658 





.68 


3 


.164 


± 





.017 


3 


.650 


± 





.025 


FLWO 


4Shooter 


3173. 


.659 


4 


.56 


3 


.210 


± 





.007 


3 


.654 


± 





.008 


MDM 


Tcmpleton 


3174. 


.664 





.71 


3 


.173 


± 





.016 


3 


.591 


± 





.023 


FLWO 


4Shooter 


3176. 


.666 


1 


.00 


3 


.172 


± 





.016 


3 


.596 


± 





.022 


FLWO 


4Shooter 


3177. 


.675 


1 


.14 


3 


.167 


± 





.018 


3 


.661 


± 





.027 


FLWO 


4Shooter 


3178. 


.675 





.63 


3 


.197 


± 





.027 


3 


.652 


± 





.040 


FLWO 


4Shooter 


3182. 


.654 





.67 


3 


.134 


± 





.043 


3 


.604 


± 





.067 


FLWO 


4 Shooter 


3183 


.653 





.66 


3 


.186 


± 





.040 


3 


.578 


± 





.057 


FLWO 


4Shooter 


3184. 


.653 





.64 


3 


.237 


± 





.036 


3 


.651 


± 





.052 


FLWO 


4Shooter 


3185. 


.652 





.60 


3 


.178 


± 





.034 


3 


.619 


± 





.050 


FLWO 


4Shooter 


3186. 


.651 





.57 


3 


.198 


± 





.039 


3 


.529 


± 





.052 


FLWO 


4Shooter 


3187. 


.651 





.56 


3 


.217 


± 





.039 


3 


.582 


± 





.054 


FLWO 


4Shooter 


3188. 


.653 





.70 


3 


.243 


± 





.057 


3 


.483 


± 





.071 


FLWO 


4Shooter 


3189. 


.656 





.66 


3 


.262 


± 





.042 


3 


.666 


± 





.061 


FLWO 


4Shooter 


3191. 


.652 





.72 


3 


.156 


± 





.042 


3 


.608 


± 





.063 


FLWO 


4Shooter 


3192. 


.660 





.81 


3 


.262 


± 





.020 


3 


.613 


± 





.027 


FLWO 


4Shooter 


3193 


.650 





.78 


3 


.254 


± 





.034 


3 


.618 


± 





.047 


FLWO 


4Shooter 


3194. 


.654 





.95 


3 


.253 


± 





.027 


3 


.571 


± 





.036 


FLWO 


4Shooter 


3195. 


.652 





.63 


3 


.208 


± 





.038 


3 


.610 


± 





.054 


FLWO 


4Shooter 


3310 


.011 





.35 


3 


.245 


± 





,041 


3 


.544 


± 





.052 


FLWO 


Minicam 


3314. 


.010 





.27 


3 


.169 


± 





.052 


3 


.576 


± 





.075 


FLWO 


Minicam 


3315 


.021 





.17 


3 


.178 


± 





.119 


3 


.782 


± 





.199 


FLWO 


Minicam 


3318. 


.964 





.80 


3 


.214 


± 





.030 


3 


.537 


± 





.039 


FLWO 


Minicam 


3321 


.023 





.55 


3 


.241 


± 





.027 


3 


.607 


± 





.036 


FLWO 


Minicam 


3335 


.040 





.60 


3 


.207 


± 





.034 


3 


.656 


± 





.049 


FLWO 


Minicam 


3336 


.026 





.48 


3 


.208 


± 





.032 


3 


.592 


± 





.044 


FLWO 


Minicam 


3341. 


.993 





.38 


3 


.230 


± 





.055 


3 


.670 


± 





.080 


FLWO 


Minicam 


3349 


.604 


4 


.16 


3 


.263 


± 





.065 


3 


.669 


± 





.093 


Wise 


TAVAS 


3350 


.012 





.42 


3 


.196 


± 





.029 


3 


.619 


± 





.040 


FLWO 


Minicam 


3350 


.998 





.72 


3 


.187 


± 





.023 


3 


.671 


± 





.033 


FLWO 


Minicam 


3352. 


.044 





.55 


3 


.212 


± 





.023 


3 


.633 


± 





.032 


FLWO 


Minicam 


3352 


.572 


2 


.39 


3 


.297 


± 





.059 


3 


.776 


± 





.092 


Wise 


TAVAS 


3353 


.038 





.90 


3 


.217 


± 





.026 


3 


.717 


± 





.038 


FLWO 


Minicam 


3354 


.033 





.79 


3 


.218 


± 





.023 


3 


.664 


± 





.031 


FLWO 


Minicam 


3354 


.920 





.89 


3 


.200 


± 





.021 


3 


.661 


± 





.030 


FLWO 


Minicam 


3359 


.000 


1 


.64 


3 


.261 


± 





.047 


3 


.876 


± 





.079 


Wise 


Tektronix 


3359. 


.920 


1 


.24 


3 


.239 


± 





.018 


3 


.667 


± 





.023 


FLWO 


Minicam 
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3360 


.000 


0. 


.40 


3 


.302 


± 


0. 


,093 


3 


.631 


± 


0. 


.124 


Wise 


Tektronix 


3377 


.879 


0. 


,71 


3 


.245 


± 


0. 


.026 


3 


.595 


± 


0. 


.034 


FLWO 


Minicam 


3378 


.945 


0. 


,43 


3 


.263 


± 


0, 


,027 


3 


.630 


± 


0. 


.036 


FLWO 


Minicam 


3379 


.878 


0. 


.42 


3 


.260 


± 


0. 


.025 


3 


.688 


± 


0. 


.035 


FLWO 


Minicam 


3380 


.000 


1. 


,33 


3 


.321 


± 


0. 


.050 


3 


.669 


± 


0. 


.067 


Wise 


Tektronix 


3380 


.606 


2. 


,78 


3 


.415 


± 





.049 


3 


.854 


± 


0. 


.075 


Wise 


TAVAS 


3381 


.581 


1 


.62 


3 


.463 


± 


0. 


.154 


3 


.993 


± 


0. 


.238 


Wise 


TAVAS 


3384 


.554 


2. 


.65 


3 


.312 


± 


0. 


.073 


3 


.680 


± 


0. 


101 


Wise 


TAVAS 


3385 


.816 


0. 


.45 


3 


.370 


± 


0. 


.014 


3 


.690 


± 


0. 


,015 


APO 


SPICam 


3385 


.951 


0. 


.63 


3 


.300 


± 


0. 


.024 


3 


.663 


± 


0. 


,031 


FLWO 


Minicam 


3386 


.453 


2. 


.74 


3 


.303 


± 


0. 


.079 


3 


.653 


± 


0, 


,114 


Wise 


TAVAS 


3387 


.527 


2. 


.92 


3 


.321 


± 


0. 


.087 


3 


.609 


± 


0. 


,109 


Wise 


TAVAS 


3387 


.960 


0. 


,69 


3 


.328 


± 


0. 


.027 


3 


.626 


± 


0. 


.033 


FLWO 


Minicam 


3399 


.018 


0. 


,30 


3 


.341 


± 


0, 


.067 


3 


.752 


± 


0. 


.096 


FLWO 


Minicam 


3399 


.714 


0. 


,25 


3 


.342 


± 


0. 


.049 


3 


.671 


± 


0. 


.064 


FLWO 


Minicam 


3402 


.990 


0. 


,27 


3 


.397 


± 


0, 


,041 


3 


.663 


± 


0. 


.052 


FLWO 


Minicam 


3403 


.412 


1. 


.99 


3 


.317 


± 


0. 


.063 


3 


.826 


± 


0. 


.101 


Wise 


TAVAS 


3406 


.322 


1, 


.32 


3 


.330 


± 


0. 


.099 


3 


.742 


± 


0, 


,146 


Wise 


TAVAS 


3408 


.375 


1. 


.96 


3 


.408 


± 


0. 


.067 


3 


.688 


± 


0, 


,086 


Wise 


m A "V T A <~"\ 

TAVAS 


3410 


.806 


0. 


.78 


3 


.330 


± 


0. 


.023 


3 


.664 


± 


0, 


.029 


FLWO 


Minicam 


3411 


.452 


2. 


.81 


(2 


.988 


± 


0. 


.065) 


3 


.489 


± 


0, 


107 


Wise 


TAVAS 


3412 


.404 


2. 


.20 


3 


.586 


± 


0. 


.092 


3 


.777 


± 


0, 


,108 


Wise 


TAVAS 


3416 


.000 


0. 


.63 


3 


.531 


± 


0. 


.103 


3 


.639 


± 


0. 


.113 


Wise 


Tektronix 


3417 


.000 


0. 


.90 


3 


.363 


± 


0. 


,072 


3 


.648 


± 


0. 


.093 


Wise 


Tektronix 


3419 


.000 


0. 


,51 


3 


.439 


± 


0. 


,073 


3 


.647 


± 


0. 


.088 


Wise 


Tektronix 


3430 


.930 


0. 


,28 


3 


.393 


± 


0, 


,043 


3 


.611 


± 


0. 


.051 


FLWO 


Minicam 


3431 


.738 


0. 


.59 


3 


.362 


± 


0. 


.028 


3 


.683 


± 


0. 


,036 


FLWO 


Minicam 


3432 


.721 


0. 


.55 


3 


.385 


± 


0. 


,028 


3 


.639 


± 


0. 


.035 


FLWO 


Minicam 


3433 


.740 


0, 


.93 


3 


.364 


± 


0. 


.025 


3 


.665 


± 


0, 


,031 


FLWO 


Minicam 


3439 


.643 


0. 


.39 


3 


.377 


± 


0. 


.013 


3 


.643 


± 


0, 


,014 


APO 


SPICam 


3439 


.719 


0. 


.71 


3 


.339 


± 


0. 


.025 


3 


.650 


± 


0. 


,032 


FLWO 


Minicam 


3441 


.615 


0. 


.34 


3 


.380 


± 


0, 


.014 


3 


.639 


± 


0, 


,015 


APO 


SPICam 


3441 


.797 


0. 


.60 


3 


.353 


± 


0. 


.030 


3 


.674 


± 


0, 


,038 


FLWO 


Minicam 


3442 


.734 


0. 


.84 


3 


.374 


± 





,026 


3 


.664 


± 


0. 


,032 


FLWO 


Minicam 


3443 


.748 


1. 


,17 


3 


.346 


± 


0. 


,029 


3 


.681 


± 


0. 


.038 


FLWO 


Minicam 


3459 


.333 


1. 


,39 


3 


.322 


± 


0. 


.068 


3 


.780 


± 


0. 


,101 


Wise 


TAVAS 


3462 


.836 





,31 


3 


.317 


± 


0. 


,037 


3 


.739 


± 


0. 


.052 


FLWO 


Minicam 


3462 


.859 


1. 


.17 


3 


.376 


± 


0. 


,013 


3 


.669 


± 


0. 


,015 


APO 


SPICam 


3463 


.722 


0. 


.55 


3 


.384 


± 


0. 


.031 


3 


.653 


± 


0. 


,038 


FLWO 


Minicam 


3464 


.719 


0. 


.81 


3 


.382 


± 


0. 


.023 


3 


.684 


± 


0, 


,029 


FLWO 


Minicam 


3466 


.251 


1, 


.65 


3 


.515 


± 


0, 


.083 


3 


.917 


± 


0. 


,119 


Wise 


TAVAS 
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Table 1 — Continued 



HJD x 2 /Ndof Image A Image B Observatory Detector 



3467 


320 


2 


26 


3 


494 


± 





063 


3 


793 


± 





083 


Wise 


TAVAS 


3468 


291 


1 


45 


3 


318 


± 





086 


3 


774 


± 





130 


Wise 


m A T T A c~\ 

TAVAS 


3469 


756 





34 


3 


396 


± 





055 


3 


766 


± 





074 


FLWO 


Minicam 


3470 


713 





47 


3 


364 


± 





026 


3 


707 


± 





033 


FLWO 


Minicam 


3471 


000 





50 


3 


457 


± 





085 


3 


746 


± 





109 


Wise 


Tektronix 


3471 


763 





32 


3 


386 


± 





029 


3 


763 


± 





039 


FLWO 


Minicam 


3472 


766 





85 


3 


397 


± 





025 


3 


709 


± 





031 


FLWO 


Minicam 


3473 


744 


1 


38 


3 


393 


± 





020 


3 


694 


± 





024 


FLWO 


Minicam 


3474 


743 





84 


3 


398 


± 





022 


3 


719 


± 





028 


FLWO 


Minicam 


3476 


739 





42 


3 


385 


± 





043 


3 


693 


± 





056 


FLWO 


Minicam 


3477 


735 





57 


3 


401 


± 





045 


3 


666 


± 





056 


FLWO 


Minicam 


3478 


690 





39 


3 


371 


± 





051 


3 


768 


± 





071 


FLWO 


Minicam 


3485 


738 





33 


3 


392 


± 





054 


3 


772 


± 





074 


FLWO 


Minicam 


3486 


707 





30 


3 


439 


± 





052 


3 


709 


± 





066 


FLWO 


Minicam 


3487 


675 


1 


08 


3 


397 


± 





025 


3 


711 


± 





032 


FLWO 


Minicam 


3489 


000 





87 


3 


384 


± 





072 


4 


182 


± 





145 


Wise 


Tektronix 


3492 


710 





46 


3 


425 


± 





038 


3 


725 


± 





049 


FLWO 


Minicam 


3494 


722 





78 


3 


412 


± 





025 


3 


771 


± 





033 


FLWO 


Minicam 


3495 


722 





28 


3 


349 


± 





059 


3 


697 


± 





080 


FLWO 


Minicam 


3496 


694 





42 


3 


415 


± 





028 


3 


704 


± 





035 


FLWO 


Minicam 


3497 


720 





43 


3 


439 


± 





031 


3 


727 


± 





039 


FLWO 


Minicam 


3498 


691 





58 


3 


422 


± 





026 


3 


732 


± 





033 


FLWO 


Minicam 


3499 


695 





73 


3 


402 


± 





025 


3 


742 


± 





032 


FLWO 


Minicam 


3500 


244 


1 


14 


3 


483 


± 





108 


3 


678 


± 





125 


Wise 


m A "X T A c\ 

TAVAS 


3501 


245 


1 


28 


3 


415 


± 





082 


3 


630 


± 





100 


Wise 


m A "y 7" A /~i 

TAVAS 


3501 


708 





66 


3 


443 


± 





044 


3 


670 


± 





053 


FLWO 


Minicam 


3501 


817 





21 


3 


458 


± 





020 


3 


735 


± 





024 


APO 


SPICam 


3502 


701 





37 


3 


411 


± 





033 


3 


746 


± 





043 


FLWO 


Minicam 


3503 


666 





74 


3 


444 


± 





028 


3 


734 


± 





035 


FLWO 


Minicam 


3504 


687 





24 


3 


515 


± 





068 


3 


702 


± 





073 


FLWO 


Minicam 


3509 


698 





31 


3 


459 


± 





054 


3 


720 


± 





068 


FLWO 


Minicam 


3510 


674 





60 


3 


472 


± 





052 


3 


715 


± 





064 


FLWO 


Minicam 


3516 


680 





53 


3 


471 


± 





031 


3 


780 


± 





040 


FLWO 


Minicam 


3520 


686 





65 


3 


562 


± 





035 


3 


778 


± 





042 


FLWO 


Minicam 


3521 


691 





51 


3 


515 


± 





034 


3 


825 


± 





044 


FLWO 


Minicam 


3522 


646 





39 


3 


498 


± 





047 


3 


746 


± 





057 


FLWO 


Minicam 


3524 


667 





63 


3 


527 


± 





034 


3 


833 


± 





043 


FLWO 


Minicam 


3525 


673 





45 


3 


518 


± 





034 


3 


854 


± 





044 


FLWO 


Minicam 


3527 


651 





70 


3 


507 


± 





034 


3 


898 


± 





047 


FLWO 


Minicam 


3530 


650 





40 


3 


569 


± 





039 


3 


893 


± 





051 


FLWO 


Minicam 


3536 


660 





29 


3 


547 


± 





072 


3 


883 


± 





096 


FLWO 


Minicam 
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Table 1 — Continued 



HJD X 2 /Ndof Image A Image B Observatory Detector 



3537 


.661 


0. 


.31 


3 


.545 


± 





.064 


3. 


.971 


± 





.092 


FLWO 


Minicam 


3538 


.652 


0. 


.43 


3 


.514 


± 





,060 


3 


.827 


± 





.078 


FLWO 


Minicam 


3541 


.654 


0. 


.28 


3 


.425 


± 





.117 


3. 


.984 


± 





.189 


FLWO 


Minicam 


3547 


.650 


0. 


.30 


3 


.595 


± 





.062 


4. 


.078 


± 


0. 


.094 


FLWO 


Minicam 


3549. 


.651 


0. 


,30 


3 


.646 


± 





.084 


4. 


.040 


± 


0. 


.118 


FLWO 


Minicam 


3655 


.027 


3 


,46 


3 


.548 


± 





.017 


3 


,983 


± 





.023 


MDM 


RETROCAM 


3656. 


.000 


5. 


,00 


3 


.621 


± 





.012 


3. 


,995 


± 





.012 


TV f"TV "ft IT 

MDM 


RETROCAM 


3664. 


.001 


1. 


,09 


3 


.553 


± 





.047 


4. 


,046 


± 





.072 


FLWO 


Keplercam 


3664 


.001 


4, 


,27 


3 


.629 


± 





.012 


4. 


,017 


± 





.014 


MDM 


RETROCAM 


3667. 


.967 


0. 


,59 


3 


.643 


± 





.049 


3. 


,967 


± 





.065 


FLWO 


Keplercam 


3668 


.006 


3. 


,51 


3 


.643 


± 





.017 


3. 


,979 


± 





.021 


"TV T T ~\ "ft I 

MDM 


RETROCAM 


3673 


.899 


0. 


,49 


3 


.635 


± 





.015 


3. 


,938 


± 





.017 


APO 


SPICam 


3674. 


.990 


0. 


,61 


3 


.600 


± 





.038 


3. 


,942 


± 





.050 


FLWO 


Keplercam 


3676 


.001 


1. 


,02 


3 


.563 


± 





.025 


3 


,928 


± 





.033 


FLWO 


Keplercam 


3676. 


.969 


0. 


,29 


3 


.473 


± 





.088 


3. 


,710 


± 


0. 


.108 


FLWO 


Keplercam 


3677. 


.013 


2. 


,75 


3 


.605 


± 





.015 


3 


.919 


± 





.018 


"ft T T "\ "ft I 

MDM 


RETROCAM 


3677 


.995 


9. 


,20 


3 


.624 


± 





.011 


3 


,911 


± 





.012 


"ft T T "\ "ft I 

MDM 


RETROCAM 


3678 


.997 


2. 


.98 


3 


.635 


± 





.012 


3. 


.914 


± 





.012 


"ft itt"\ ft it 

MDM 


n t~~\ m n /~\ /"I A ~\ it 

RETROCAM 


3679 


.013 


1. 


,03 


3 


.583 


± 





.034 


3. 


.881 


± 





.043 


FLWO 


Keplercam 


3679. 


.995 


0. 


.60 


3 


.581 


± 





,038 


3. 


.947 


± 





.051 


FLWO 


Keplercam 


3681. 


.004 


1, 


,26 


3 


.553 


± 





.024 


3. 


.949 


± 





.032 


FLWO 


Keplercam 


3684 


.012 


0. 


,55 


3 


.650 


± 





.052 


3, 


,856 


± 





.061 


FLWO 


Keplercam 


3685 


.017 


0. 


,70 


3 


.577 


± 





.033 


3. 


.944 


± 





.045 


FLWO 


Keplercam 


3686 


.564 


4. 


.61 


3 


.476 


± 





.066 


3 


,829 


± 





.094 


Wise 


m A "X t A n 

TAVAS 


3686 


.892 


1. 


.46 


3 


.550 


± 





.026 


3 


.891 


± 





.035 


Palomar 


SITc 


3687 


.024 


0. 


.89 


3 


.551 


± 





.035 


3 


.931 


± 





.048 


FLWO 


Keplercam 


3688. 


.016 


1. 


.28 


3 


.589 


± 





.029 


3 


.877 


± 





.037 


b LWO 


TV 1 

Keplercam 


3688. 


.929 


1. 


.86 


3 


.602 


± 





.017 


3 


.889 


± 





.020 


MDM 


RETROCAM 


3689 


.034 


0. 


,54 


3 


.575 


± 





.038 


3 


.862 


± 





.048 


FLWO 


Keplercam 


3690 


.002 


0. 


.43 


3 


.510 


± 





.067 


3. 


,904 


± 





.095 


FLWO 


Keplercam 


3691 


.019 


0. 


.56 


3 


.532 


± 





.054 


3, 


,767 


± 





.067 


FLWO 


Keplercam 


3691 


.898 


1. 


.70 


3 


.570 


± 





.031 


3. 


,929 


± 





.042 


Palomar 


SITe 


3692 


.018 


0. 


,40 


3 


.607 


± 





.056 


3. 


,869 


± 





.070 


FLWO 


Keplercam 


3693 


.022 


0. 


,38 


3 


.535 


± 





.054 


3 


,949 


± 





.077 


FLWO 


Keplercam 


3693 


.865 


0. 


,59 


3 


.588 


± 





.049 


3 


,856 


± 





.062 


Palomar 


SITe 


3693 


.927 


0. 


,28 


3 


.607 


± 





.023 


3 


,894 


± 





.028 


APO 


SPICam 


3694 


.012 


0. 


.61 


3 


.557 


± 





.048 


3 


,986 


± 





.069 


FLWO 


Keplercam 


3694 


.862 


0. 


,77 


3 


.654 


± 





.046 


3 


.977 


± 





.061 


Palomar 


SITe 


3698 


.919 


3. 


.60 


3 


.640 


± 





.012 


3. 


.945 


± 





.013 


MDM 


RETROCAM 


3700 


.601 


0. 


.69 


-0.13 


± 


1. 


048 


(-0.51 


± 


0. 


926) 


Wise 


TAVAS 


3700 


.927 


1. 


.32 


3.592 


± 





.021 


3.915 


± 





.027 


FLWO 


Keplercam 
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Table 1 — Continued 



HJD x 2 /N do f Image A Image B Observatory Detector 



3700 


.997 


1.99 


3 


.618 


± 





.012 


(3. 


.851 


± 





.013) 


MDM 


RETROCAM 


3701 


.581 


0.87 


3, 


,002 


± 





.238 


(3. 


.054 


± 





.238) 


Wise 


m A T T A c~\ 

TAVAS 


3702 


.560 


2.57 


3. 


,456 


± 





.104 


3 


.595 


± 





.115 


Wise 


m A t r A i~i 

TAVAS 


3706 


.012 


3.42 


3 


,611 


± 





.019 


3 


.934 


± 





.024 


FLWO 


Kcplercam 


3707 


.581 


3.31 


3, 


,480 


± 





.067 


4. 


,099 


± 





.120 


Wise 


m A T 7" A C\ 

TAVAS 


3708 


.849 


0.45 


3 


,539 


± 





.045 


3 


,992 


± 





.070 


Palomar 


SITc 


3708 


.984 


0.99 


3, 


,629 


± 





.029 


3. 


,927 


± 





.036 


FLWO 


Kcplercam 


3709 


.996 


0.79 


3, 


.561 


± 





.030 


3. 


.944 


± 





.041 


FLWO 


Keplercam 


3710 


.561 


0.94 


3 


.558 


± 





.186 


4, 


.183 


± 





.311 


Wise 


m A T T A r~1 

TAVAS 


3710 


.896 


0.56 


3, 


.613 


± 





.013 


3 


.934 


± 





.015 


APO 


SPICam 


3710 


.971 


1.46 


3 


.618 


± 





.026 


3 


.964 


± 





.034 


FLWO 


Keplercam 


3711 


.943 


8.99 


3. 


,602 


± 





.013 


3 


.979 


± 





.014 


FLWO 


Kcplercam 


3712 


.545 


5.05 


3 


,592 


± 





.063 


3 


,917 


± 





.085 


Wise 


m A t 7 A n 

TAVAS 


3712 


.871 


1.21 


3. 


,533 


± 





.022 


4. 


.016 


± 





.032 


FLWO 


Kcplercam 


3712 


.891 


4.05 


(3. 


.667 


± 





.012) 


4, 


,005 


± 





.012 


MDM 


RETROCAM 


3713 


.483 


3.09 


3, 


,635 


± 





.059 


4. 


.142 


± 





.094 


Wise 


TAVAS 


3714 


.951 


0.54 


3 


603 


± 





.044 


4. 


.008 


± 





.062 


FLWO 


Kcplercam 


3718 


.765 


0.56 


3 


,581 


± 





.054 


4, 


.023 


± 





.079 


Palomar 


SITc 


3725 


.965 


5.92 


3 


.576 


± 





.016 


4. 


.032 


± 





.021 


FLWO 


Kcplercam 


3726 


.932 


5.41 


3, 


,588 


± 





.014 


4, 


.095 


± 





.019 


FLWO 


Keplercam 


3728 


.935 


0.62 


3, 


.593 


± 





.039 


4. 


.069 


± 





.059 


FLWO 


Keplercam 


3729 


.951 


17.11 


3 


,609 


± 





.011 


4. 


.090 


± 





.012 


TV IT 1 "\ TV ,T 

MDM 


RETROCAM 


3729 


.984 


1.96 


3 


.581 


± 





.021 


4. 


,059 


± 





.030 


FLWO 


Kcplercam 


3730 


.951 


0.65 


3 


,549 


± 





.049 


4, 


.050 


± 





.076 


FLWO 


Kcplercam 


3731 


.854 


0.60 


3. 


.556 


± 





.035 


4. 


.038 


± 





.053 


FLWO 


Kcplercam 


3732 


.464 


3.27 


3 


.502 


± 





.053 


4. 


.231 


± 





.105 


Wise 


TAVAS 


3734 


.457 


3.01 


3. 


,594 


± 





.067 


3 


.941 


± 





.092 


Wise 


TAVAS 


3734 


.918 


12.36 


3 


,633 


± 





.011 


4. 


.114 


± 


0. 


.012 


MDM 


RETROCAM 


3736 


.018 


3.92 


3, 


.590 


± 





.011 


4. 


.056 


± 





.012 


MDM 


RETROCAM 


3736 


.902 


3.77 


3, 


.629 


± 





.011 


4. 


.073 


± 





.012 


MDM 


RETROCAM 


3737 


.923 


3.42 


3, 


.621 


± 





.012 


4. 


.087 


± 





.014 


MDM 


RETROCAM 


3738 


.961 


2.73 


3, 


.603 


± 





.014 


4, 


.082 


± 





.018 


MDM 


RETROCAM 


3740 


.032 


0.60 


3 


.540 


± 





.030 


4. 


.047 


± 





.045 


FLWO 


Keplercam 


3740 


.910 


3.10 


3. 


661 


± 





.011 


4, 


.121 


± 





.012 


MDM 


RETROCAM 


3741 


.012 


0.61 


3, 


,488 


± 





.056 


3 


.982 


± 





.088 


FLWO 


Keplercam 


3741 


.885 


1.76 


(3, 


680 


± 





.011) 


(4- 


.152 


± 





.012) 


MDM 


RETROCAM 


3742 


.992 


1.50 


3 


,565 


± 





.019 


4. 


.109 


± 





.028 


FLWO 


Kcplercam 


3743 


.772 


5.24 


3. 


,519 


± 





.018 


4, 


.087 


± 





.028 


Palomar 


SITe 


3743 


.976 


19.51 


3, 


.660 


± 





.011 


4. 


.119 


± 





.012 


MDM 


RETROCAM 


3744 


.962 


379.63 


3, 


.626 


± 





.052 


4. 


.081 


± 





.052 


MDM 


RETROCAM 


3745 


.931 


17.41 


3, 


.600 


± 





.014 


4. 


.088 


± 





.017 


MDM 


RETROCAM 
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3745 


.970 


0.70 


3 


.619 


± 





.032 


4, 


.054 


± 





.046 


FLWO 


Keplercam 


3747 


.024 


8.47 


3. 


,653 


± 





.012 


4. 


,088 


± 


0. 


.012 


T\ IT 1 \ TV ,T 

MDM 


RETROCAM 


3749 


.000 


2.14 


3. 


.651 


± 





.010 


4, 


.074 


± 





.014 


TV IT 1 \ TV ,T 

MDM 


Echelle 


3752 


.918 


1.12 


3 


643 


± 





.017 


4. 


,074 


± 





.024 


MDM 


Echclle 


3753 


.950 


1.71 


3 


,610 


± 





.020 


4. 


,000 


± 





.027 


FLWO 


Keplercam 


3755 


.878 


2.61 


3. 


,650 


± 





.012 


(4. 


.162 


± 


0. 


.013) 


TV I 1 "\ TV ,T 

MDM 


T" V 1 1 m T" V / ~\ S * A TV T 

RETROCAM 


3755 


.897 


0.41 


3 


.591 


± 





.041 


4, 


.044 


± 





.062 


FLWO 


Keplercam 


3756 


.929 


3.27 


3 


.654 


± 





.012 


4, 


.057 


± 





.012 


T\ 1 1 "\ TV T" 

MDM 


RETROCAM 


3757 


.884 


14.89 


3 


.657 


± 





.011 


4, 


.047 


± 





.011 


TV I 1 "\ TV ,T 

MDM 


T" V 1 1 m T" V S \ S \ A TV T 

RETROCAM 


3757 


.892 


0.74 


3 


.605 


± 





.034 


4, 


.035 


± 





.048 


FLWO 


Keplercam 


3758 


.916 


10.35 


3 


.657 


± 





.011 


4, 


.034 


± 





.011 


MDM 


RETROCAM 


3758 


.937 


0.55 


3. 


,618 


± 





.032 


4, 


.016 


± 





.043 


FLWO 


Keplercam 


3761 


.985 


0.56 


3 


,648 


± 





.095 


3 


,989 


± 





.128 


Palomar 


SITe 


3762 


.389 


2.23 


3. 


,575 


± 





.086 


4. 


.034 


± 





.130 


Wise 


TAVAS 


3764 


.885 


0.50 


3. 


,646 


± 





.043 


4, 


.056 


± 





.060 


FLWO 


Keplercam 


3766 


.051 


0.58 


3 


665 


± 





.034 


4, 


.084 


± 





.048 


FLWO 


Keplercam 


3766 


.906 


1.01 


3. 


,685 


± 





.026 


4, 


.033 


± 





.034 


FLWO 


Keplercam 


3767 


.858 


5.20 


3 


,671 


± 





.013 


4. 


.065 


± 





.015 


FLWO 


Keplercam 


3768 


.842 


1.03 


3 


,693 


± 





.030 


3. 


,971 


± 





.037 


FLWO 


Keplercam 


3769 


.909 


1.19 


3 


,616 


± 





.054 


4. 


,056 


± 





.079 


FLWO 


Keplercam 


3770 


.792 


2.10 


3, 


.700 


± 





.011 


4. 


,038 


± 





.013 


MDM 


8K 


3770 


.908 


0.86 


3, 


,666 


± 





.027 


4. 


,078 


± 





.038 


FLWO 


Keplercam 


3771 


.393 


3.88 


3 


.830 


± 





.082 


4, 


.045 


± 





.099 


Wise 


m A t t A 

TAVAS 


3771 


.772 


5.09 


3. 


,642 


± 





.016 


3 


.940 


± 





.020 


MDM 


8K 


3771 


.870 


0.63 


3 


,748 


± 





.046 


4. 


.074 


± 





.061 


FLWO 


Keplercam 


3772 


.435 


1.70 


3 


612 


± 





.083 


3 


.915 


± 





.108 


Wise 


TAVAS 


3772 


.758 


0.83 


3. 


,633 


± 





.030 


4. 


.150 


± 





.047 


MDM 


8K 


3772 


.922 


0.57 


3. 


,754 


± 





.041 


4. 


.150 


± 





.058 


FLWO 


Keplercam 


3773 


.722 


1.70 


3. 


,699 


± 





.014 


4. 


.044 


± 





.018 


MDM 


8K 


3775 


.952 


0.54 


3. 


.729 


± 





.035 


3. 


959 


± 





.042 


Palomar 


SITe 


3776 


.938 


5.19 


3. 


.691 


± 





.022 


3, 


,997 


± 





.028 


MDM 


8K 


3786 


.858 


0.49 


3 


.699 


± 





.039 


4. 


065 


± 





.053 


FLWO 


Keplercam 


3787 


.802 


0.35 


3. 


.721 


± 





.089 


4. 


,069 


± 





.120 


FLWO 


Keplercam 


3788 


.859 


1.36 


3. 


.676 


± 





.022 


4. 


.084 


± 





.029 


FLWO 


Keplercam 


3790 


.773 


1.79 


3. 


,704 


± 





.043 


4. 


.001 


± 





.055 


FLWO 


Keplercam 


3791 


.791 


1.15 


3. 


,648 


± 





.029 


4. 


.039 


± 





.039 


FLWO 


Keplercam 


3792 


.946 


0.39 


3. 


,730 


± 





.112 


4. 


.290 


± 





.180 


Palomar 


SITe 


3793 


.794 


1.21 


3. 


.675 


± 





.028 


4. 


.016 


± 





.037 


FLWO 


Keplercam 


3794 


.750 


0.64 


3 


.681 


± 





.038 


4. 


.071 


± 





.052 


FLWO 


Keplercam 


3795 


.306 


4.37 


3 


,668 


± 





.069 


4. 


,261 


± 





.118 


Wise 


TAVAS 


3797 


.286 


2.39 


3. 


.643 


± 





.068 


4. 


,005 


± 





.092 


Wise 


TAVAS 
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3797. 


.797 





,47 


3. 


,599 


± 





.047 


4. 


,066 


± 





.069 


FLWO 


Keplercam 


3798. 


.375 


4 


.00 


3. 


,576 


± 





.043 


3. 


,928 


± 





.059 


Wise 


m A t t A <~1 

TAVAS 


3798. 


.841 


1 


.42 


3. 


.691 


± 





.031 


3. 


.951 


± 





.038 


FLWO 


Keplercam 


3799 


.364 


5 


.65 


3. 


696 


± 





.040 


4. 


,009 


± 





.053 


Wise 


m A "Y 7" A /~1 

TAVAS 


3799 


.854 





.73 


3. 


,682 


± 





.034 


3. 


.972 


± 





.042 


FLWO 


Keplercam 


3800. 


.037 





,69 


3. 


,640 


± 





.047 


3 


,883 


± 





.058 


Palomar 


SITe 


3800. 


.791 





.29 


3, 


,602 


± 





.091 


4. 


,060 


± 





.135 


FLWO 


Keplercam 


3815 


,828 





.71 


3, 


,712 


± 





.043 


4, 


.013 


± 





.055 


FLWO 


Keplercam 


3816. 


330 


2 


.44 


3, 


,675 


± 





.050 


4. 


.067 


± 





.070 


Wise 


TAVAS 


3817. 


,854 





.51 


3, 


,698 


± 





.051 


3. 


.985 


± 





.067 


FLWO 


Keplercam 


3817. 


,933 





.24 


3, 


.818 


± 





.146 


4, 


.035 


± 





.175 


Palomar 


SITe 


3818 


,272 


4 


.80 


3 


,709 


± 





.111 


3 


.759 


± 





.117 


Wise 


m A t t A <~1 

TAVAS 


3818. 


,843 





.86 


3, 


,667 


± 





.048 


4. 


.002 


± 





.064 


FLWO 


Keplercam 


3819. 


,721 


1 


.05 


3. 


.815 


± 





.059 


3. 


.925 


± 





.065 


FLWO 


Keplercam 


3820. 


,240 


6 


.87 


3, 


,637 


± 





.062 


4. 


.041 


± 





.090 


Wise 


TAVAS 


3820. 


,720 





.59 


3, 


,682 


± 





.045 


3. 


.958 


± 





.056 


FLWO 


Keplercam 


3822. 


.247 


2 


.33 


3. 


,671 


± 





.082 


4. 


.052 


± 





.114 


Wise 


TAVAS 


3823. 


696 


1 


.49 


3, 


,698 


± 





.027 


3. 


.980 


± 





.034 


FLWO 


Keplercam 


3827. 


,816 





.57 


3, 


,722 


± 





.037 


3. 


.940 


± 





.045 


FLWO 


Keplercam 


3827. 


,919 





.66 


3, 


,751 


± 





.028 


3. 


.937 


± 





.032 


Palomar 


SITe 


3832. 


,752 





.55 


3, 


.674 


± 





.073 


4, 


,092 


± 





.106 


FLWO 


Keplercam 


3837. 


,715 





.47 


3, 


,712 


± 





.070 


3, 


.937 


± 





.085 


FLWO 


Keplercam 


3841 


.669 





.67 


3. 


,747 


± 





,042 


3. 


,923 


± 





.049 


FLWO 


Keplercam 


3843. 


.674 


1 


.23 


3. 


,654 


± 





.034 


3. 


,905 


± 





.041 


FLWO 


Keplercam 


3844. 


.642 





.81 


3. 


,715 


± 





.035 


3. 


,952 


± 





.042 


FLWO 


Keplercam 


3845. 


.668 





.75 


3. 


,668 


± 





.035 


3. 


,864 


± 





.041 


FLWO 


Keplercam 


3846. 


.265 


1 


.84 


3. 


,703 


± 





.063 


3. 


,841 


± 





.071 


Wise 


TAVAS 


3846, 


.644 





.94 


3. 


,702 


± 





.033 


3 


.895 


± 





.038 


FLWO 


Keplercam 


3848. 


665 


1 


.15 


3, 


.659 


± 





.030 


3. 


.960 


± 





.038 


FLWO 


Keplercam 


3849, 


.667 


1 


.41 


3, 


.658 


± 





.037 


3. 


.864 


± 





.044 


FLWO 


Keplercam 


3850. 


696 


1 


.04 


3, 


.689 


± 





.035 


3. 


.844 


± 





.039 


FLWO 


Keplercam 


3851 


.707 





.92 


3, 


.690 


± 





.032 


3. 


.886 


± 





.037 


FLWO 


Keplercam 


3852. 


.837 





.38 


3, 


,701 


± 





.065 


3. 


.958 


± 





.082 


FLWO 


Keplercam 


3854. 


.730 





.62 


3. 


,707 


± 





.030 


3. 


.941 


± 





.036 


FLWO 


Keplercam 


3855. 


.650 


2 


.11 


3. 


,668 


± 





.025 


3 


,897 


± 


0. 


.030 


FLWO 


Keplercam 


3856. 


.695 


5 


.33 


3, 


.659 


± 





.015 


3. 


.883 


± 





.016 


FLWO 


Keplercam 


3857. 


.822 





.55 


3. 


,711 


± 





.045 


3 


,874 


± 





.052 


FLWO 


Keplercam 


3858. 


.643 


1 


.17 


3. 


,629 


± 


0. 


.029 


3. 


.889 


± 





.036 


FLWO 


Keplercam 


3859 


.641 


1 


.63 


3, 


,614 


± 





.032 


3. 


.929 


± 





.041 


FLWO 


Keplercam 


3861. 


.765 





.57 


3, 


,736 


± 





.045 


3, 


.966 


± 





.055 


FLWO 


Keplercam 


3873. 


.245 


3 


.62 


3, 


.645 


± 





.057 


4. 


.150 


± 





.089 


Wise 


TAVAS 
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3875. 


.259 





.75 


3, 


.605 


± 





.107 


3 


.992 


± 





.151 


Wise 


TAVAS 


3878, 


.251 


2 


.27 


3. 


,686 


± 





.054 


4. 


,044 


± 





.074 


Wise 


TAVAS 


3879. 


.245 


1 


.27 


3, 


.676 


± 





.068 


4. 


,071 


± 





.095 


Wise 


TAVAS 


3881. 


.706 


1 


.46 


3, 


,669 


± 





.025 


4. 


.035 


± 





.033 


FLWO 


Keplercam 


3882. 


666 


1 


.81 


3. 


.664 


± 





.024 


4. 


,021 


± 





.031 


FLWO 


Keplercam 


3883. 


.654 


1 


.35 


3, 


.608 


± 





.026 


4. 


,074 


± 





.038 


FLWO 


Keplercam 


3884. 


.724 


1 


.76 


3, 


.630 


± 





.010 


4, 


.004 


± 





.013 


MDM 


8K 


3885. 


653 


1 


.77 


3, 


613 


± 





.027 


4. 


.009 


± 





.037 


FLWO 


Keplercam 


3885. 


.720 


3 


.67 


3, 


.648 


± 





.011 


4, 


.038 


± 





.013 


MDM 


8K 


3886, 


.646 


1 


.99 


3. 


.595 


± 





.011 


4, 


.003 


± 





.014 


MDM 


8K 


3887, 


.646 


1 


.06 


3, 


.655 


± 





.014 


4, 


.010 


± 





.018 


MDM 


8K 


3890, 


.729 





.40 


3, 


,546 


± 





.058 


3 


,916 


± 





.080 


FLWO 


Keplercam 


3896 


.660 





.37 


3, 


.554 


± 





.057 


3, 


.980 


± 





.083 


FLWO 


Keplercam 


3899 


.680 





.33 


3, 


.466 


± 





.074 


3 


,941 


± 





.113 


FLWO 


Keplercam 


3903 


.708 





.57 


3. 


.629 


± 





.040 


3, 


,954 


± 





.053 


FLWO 


Keplercam 


3907, 


.679 





.84 


3 


.619 


± 





.030 


3, 


,995 


± 





.040 


FLWO 


Keplercam 



Note. — The Heliocentric Julian Days (HJD) are days relative to HJD= 2450000. 
The x 2 /Ndof column indicates how well our photometric model fit the imaging 
data. When x 2 /Ndof we rescale the photometric errors presented in this Table by 
(x 2 /N do f) x / 2 before carrying out the time delay analysis to reduce the weight of im- 
ages that were fit poorly. The image A and B magnitudes are relative to the comparison 
stars (see text). The 16 magnitudes enclosed in parentheses are not used in the time 
delay estimates. 
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